ON HIGH SPOTS OF THE FUNDAMENTAL SLOSHING 
EIGENFUNCTIONS IN AXIALLY SYMMETRIC DOMAINS 

TADEUSZ KULCZYCKI AND MATEUSZ KWASNICKI 

Abstract. We investigate the classical eigenvalue problem that arises in hydrody- 
namics and is referred to as the sloshing problem. It describes free liquid oscillations 
in a liquid container IV C I 3 . We study the case when W is an axially symmetric, 
convex, bounded domain satisfying the John condition. The Cartesian coordinates 
(x, y, z) are chosen so that the mean free surface of the liquid lies in (x, z)-plane and 
y-axis is directed upwards (y-axis is the axis of symmetry) . Our first result states that 
the fundamental eigenvalue has multiplicity 2 and for each fundamental eigenfunction 
if there is a change of x, z coordinates by a rotation around y-axis so that <p is odd in 
x- variable. 

The second result of the paper gives the following monotonicity property of the 
fundamental eigenfunction ip. If (p is odd in cc-variable then it is strictly monotonic in 
x-variable. This property has the following hydrodynamical meaning. If liquid oscil- 
lates freely with fundamental frequency according to ip then the free surface elevation 
of liquid is increasing along each line parallel to x-axis during one period of time and 
decreasing during the other half period. The proof of the second result is based on 
the method developed by D. Jerison and N. Nadirashvili for the hot spots problem for 
Neumann Laplacian. 



1. Introduction 

Linear water-wave theory is a widely-used approach that allows to determine the 
frequencies and modes of free oscillations of a liquid in a container. Such oscillations 
exist provided the liquid's upper surface is free and, in the framework of this theory, one 
obtains their frequencies and modes from a mixed Steklov problem. The latter involves 
a spectral parameter in the boundary condition on the free surface. This boundary 
value problem (usually referred to as the sloshing problem) has been the subject of a 
great number of studies over 250 years (see [10] for a historical review and [3], [13] , 
[IB] ; [22], [23], [21] for some of recent literature). It is also worth pointing out here 
that other Steklov type eigenvalue problems have attracted considerable attention in 
last years. For some of these developments, see e.g. [15], [IT], [2], [9]. 

Recently, the question of the so-called 'high spots' defined by sloshing eigenf unctions 
corresponding to the fundamental eigenvalue attracted the authors' attention. This 
question is not only similar, but closely related to the long-standing 'hot spots' conjec- 
ture of J. Rauch. (It is worth mentioning that a substantial progress has been achieved 
in studies of this conjecture for the Neumann Laplacian during the past decade; see, for 
example, the works [31], [TJ, [20], [7], [BJ.) Roughly speaking, the question about high 
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spots concerns monotonicity properties of fundamental sloshing eigenfunctions (see sub- 
section 1.3 for a detailed description). Several results about the location of high spots 
were proved in [26] and [27] . One of them deals with the fundamental eigenfunction 
(it is unique up to a non-zero factor) of the two-dimensional sloshing problem in the 
case when the domain's top interval is the one-to-one orthogonal projection of the bot- 
tom. The other one treat fundamental eigenfunctions in troughs (their cross-sections 
are subject to the same condition), and some vertical axisymmetric containers. More- 
over, it was shown in [26] that for vertical-walled containers with horizontal bottom the 
question about high spots is equivalent to the hot spots conjecture. 

The aim of this paper is to study the location of high spots for fundamental eigenfunc- 
tions satisfying the three-dimensional sloshing problem in axially symmetric domains 
of rather general shape. It occurs that the method, which can be briefly characterised 
as the method of domain's deformation (it was developed by D. Jerison and N. Nadi- 
rashvili in [20] in order to prove the hot spots conjecture for domains with two axes of 
symmetry) is adaptable for our purpose. The result demonstrating that eigenvalues and 
eigenfunctions of the sloshing problem depend continuously on the domain deformation 
are of interest in itself. 



1.1. Sloshing problem. First we formulate the three-dimensional sloshing problem in 
its general form. 

Let an inviscid, incompressible, heavy liquid occupy a three-dimensional container 
bounded from above by a free surface, which in its mean position is a simply connected 
two-dimensional domain of finite diameter. Let Cartesian coordinates (x, y, z) be chosen 
so that the mean free surface lies in the (x, ,z)-plane and the ?/-axis is directed upwards. 
The surface tension is neglected on the free surface, and we assume the liquid motion 
to be irrotational and of small amplitude. The latter assumption allows us to linearise 
boundary conditions on the free surface and this leads to the following boundary value 
problem for ip(x,y, z) — the velocity potential of the flow with a time-harmonic factor 
removed: 



Aip = in W, (1.1) 

^T = v<P on F, (1.2) 

oy 

^ = on B. (1.3) 

on 

Here W C {(x, y, z) G R 3 : y < 0} is the domain which is supposed to be a bounded Lip- 
schitz domain. The boundary dW consists of a two-dimensional domain F C {(x, 0, z) : 
x, z G R} referred to as the free surface (we assume that F is a bounded Lipschitz do- 
main) and B = dW \ F, the rigid container's bottom. In the whole paper we will refer 
to a domain W with the above geometric properties as to a liquid domain. Throughout 
the article, 4- denotes the normal derivative at dW, which is well-defined for almost 
every (with respect to the surface me asur e) point of dW. We understand that if is a 

JL 

On 



continuous function on W, and that (|l.3|) is satisfied for all points on B for which " 
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Figure 1. (a) An axisymmetric container W obtained as rotation of a 
domain D. (b) A domain D. 



is defined. We remark that in condition (1.2), the coefficient v = u 2 /g is the spec- 



tral parameter which involves the radian frequency uj of liquid's oscillations and the 
acceleration due to gravity g. 



The zero eigenvalue obviously exists for the problem (1.1 )— ( 1.3 ), but we exclude it 
with the help of the following orthogonality condition: 



0. 



;i-4) 



It has been known since the 1950s that the problem (1.1 )— ( 1.4 ) has a discrete spectrum; 
that is, there exists a sequence of eigenvalues 

< ui < v 2 < ■ ■ ■ < v n < . . . , (1.5) 



each having a finite multiplicity equal to the number of repetitions in (1.5), and such 



that v n — > oo as n — > oo. These eigenvalues can be found by means of the variational 
principle (see, for example, [30]), corresponding to the following Rayleigh quotient: 



chp\2 
dx> 



+ 



'chp\2 
> dy ' 



+ 



■ dtp 
• dz 



) 2 ] 



1.6) 



where ip is in the Sobolev space H l {W) and satisfies (1.4). Thus v\ is equal to the 
minimum of this quotient over the subspace of the Sobolev space H l {W) which consists 
of functions satisfying (1.4); the corresponding eigenfunction delivers the value v\ to the 
quotient. The eigenfunctions <p n (n = 1, 2, . . . ) belong to H^-iW) and form a complete 
system in an appropriate Hilbert space (see, for example, [22J). 

1.2. Axisymmetric containers. Now we turn to the problem of sloshing in axisym- 
metric containers. It is convenient to introduce the cylindrical coordinates (r, 6, y) so 
that 

x = rcos#, z = rsm9. (1-7) 

and to take the y-axis as the axis of symmetry for W. In this case F is typically a disc 
on xz-plane (see Figure [ija)). Moreover, we will consider W as obtained by rotation 
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of a domain D adjacent to both axes in the ry-plane. For such liquid domains, we 
write W = W(D). It is convenient to think of D as the cross-sections of W along 
the half-plane 6 = 0. By F(D) and B(D) we denote the cross-sections of F and B, 
respectively, while R(D) is the part of dD located on the y-axis (see Figure [lib)). 
It is clear that the ansatz 



ip = ip(r, y) cos(m6>), 



or 



ip = ip(r, y) sin(m#), m = 0,1,2, 



1. 



where ip is bounded near R(D), reduces the eigenvalue problem (1.1)-(1.4) in W to the 
following sequence of boundary value problems: 

,2 



d 2 ip 1 dip d 2 ip 
dr 2 r dr dy 2 



m 



ip = 



dip 
dy 
dip 
dn 



vip 



These relations must hold for all m, and by (1.4), for m 



in D, (1.9) 

on F(D), (1.10) 

on B(D). (1.11) 
we also have that 



ip(r, y)r dr dy = 0. 



(1.12) 



F(D) 



The above reduction was applied by many authors, see, in particular, [291 P- 56], 
p. 294], [131 formulas (13), (14)]. 



The variational method guarantees that for every m 



lem ( 1.9 )-( 1.12) has a sequence of eigenvalues 
< JVi 



< ^m,2 < V: 



m,3 



< 



rn 



0,1,2, 



0,1,2, 



the spectral prob- 



(1.13) 



and by ip m k, m > 0, k > 1, we denote the double sequence of corresponding eigen- 



functions. Every eigenvalue in (1.13) has a finite multiplicity equal to the number of 



repetitions; moreover, for every m > we have that v m k — > oo as k — » oo. 



It is clear that, the sequence of eigenvalues {i/ n }™ =1 of problem (1.1 )— ( 1.4 ) for W 
coincides with the double sequence {v m ,k\m=o k=i-> w hh every number v m ^ repeated 



twice when m > 1. Thus the sequence of problems ( 1.9 )-( 1.12) is equivalent to the 
original sloshing problem. 

1.3. Statement of results. We need the following definition. We say that a liquid 
domain W satisfies the John condition when W <Z F x (— oo,0). 

The following theorem is the main result of this paper (cf. [201 Theorem 1.1], [261 
Theorem 2.1], [27| Theorem 3.1]). 



Theorem 1.1. Let us consider the sloshing problem (1.1 - l.Jj.). Assume that a liquid 



domain W is an axisymmetric, convex, bounded domain, satisfying the John condition. 
We consider W as obtained by rotation of a domain D (see Figure^. Then we have: 

(i) The fundamental eigenvalue V\ equals and has multiplicity 2. Two linearly 
independent eigenf unctions corresponding to v\ are given in cylindrical coordi- 
nates by ipi t i(r, y) cos 8 and ipi,i(r, y) sin#. 
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(ii) After multiplication by ±1 we may assume that ip^i > on D. We have 

Wit ^ n ^ n n 

— — > 0, — — > on D. 

or ay 

(iii) Let <p(x, y, z) be the eigenfunction corresponding to v\ which is odd in x variable 
(in cylidrical coordinates <p equals i>i,i(r, y) cos9). Denote W + = {(x,y,z) G 
W : x > 0}. After multiplication by ±1 we may assume that > on W+. 
Then 

^ > on W, > on W+. 

ox ay 

For any other eigenfunction tp corresponding to v\ there is a change of x, z 

coordinates by a rotation around y-axis so that tp is odd in x variable. 

The third part of the theorem has the following hydrodynamical meaning. If liquid 
oscillates freely with the fundamental frequency v\ then at every moment the free- 
surface elevation of liquid is proportional to the fundamental eigenfunction ip(x, 0, z) 
(see e.g. [2H])- If we assume that ip(x,y,z) is odd in x variable then the elevation is 
increasing along each line parallel to x axis during one half-period of liquid oscillation 
and decreasing during the other half-period. In particular when the free surface is 
F = {(x, y, z) : x 2 + z 2 < r%, y = 0} then the elevation has its maximum at (r , 0, 0) 
and minimum at (— ro,0,0) during one half-period of oscillation, whereas during the 
other half-period the maximum and minimum values exchange places with one another. 
This is the reason to call this property the 'high spots' theorem. 

One could ask whether the assumption that W satisfies the John condition is neces- 



sary in Theorem 1.1 It occurs that if W does not satisfy the John condition then the 



monotonicity property of ip\ i does not necessarily hold. More precisely we have 



Proposition 1.2. Let us consider the sloshing problem (1.1 - l.Jf). Assume that W 
is an axisymmetric liquid domain for which F is a disk and B is a C 2 surface. We 
consider W as obtained by rotation of a domain D. After multiplication by ±1 we may 
assume that ipi t \ > on D. If the angle between F(D) and B(D) at the point where 
F(D) and B(D) meet (see Figure^b)) is bigger than it/2 and smaller than tt then ipi^ 
attains maximum in the interior of F(D), and changes the sign in D. 

One could also ask whether assumption about convexity of W is necessary in Theorem 
|1.1| Indeed we will show monotonicity property of ipx t i for slightly more general class 



of domains, see Definition 3.1 and Theorem 3.3 However for this class of domains we 
were not able to prove that v% = v\±. 

Although numerical results strongly suggest that Theorem |1.1| (i) should hold the 
proof of Theorem [lT] (i) is far from being trivial. The most difficult part of this proof is 
to show that v\ is not ^o,i which is the smallest eigenvalue corresponding to an axially 
symmetric eigenfunction. The proof of Theorem |1.1| (i) is based on results by Troesch 
[35] obtained by inverse methods. 

1.4. Organization of the paper. In Section 2 we prove Theorem |1.1| (i). The rest 
of the paper deals with monotonicity properties of fundamental eigenfunctions. We 
use methods from [20J, which may be briefly described as deformations of domains. In 
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Section 3 we first define a new class of domains W and formulate monotonicity properties 
of ipi^i for this class (see Theorem 3.3). Then we prove continuos dependence of u^i 
and ipi t i under certain variations of the domain. In Section 4 we prove monotonicity 
properties of ipi t i for a special class of piecewise smooth domains. In Section 5 we pass 
to the limit to obtain the same result for class W. As a conclusion we obtain Theorem 
[T7T1 (ii) and (iii). 

Throughout the article, except Definition |3.4[ Lemma 3J3 and Lemma 3^ we only 
study axisymmetric water domains W = W(D), with free surface F being a disk. We 
switch freely between Cartesian and cylindrical coordinate systems. By scaling, it will 
be often sufficient to consider the case when F is the unit disk {(x,y,z) : x 2 + y 2 < 
1, y — 0}. By a standard argument, ipi tl does not change sign in D, and v x> \ < V\^. 
With no loss of generality, we assume that ^1,1 is positive on D. We frequently use the 
continuity of ip\ y \ on D and smoothness of m D. 



2. Fundamental eigenvalue 

In this section we prove Theorem (i). The most difficult element of the proof of 
this theorem is to exclude the possibility that v\ equals z^i- 

At first we need to introduce an auxiliary Dirichlet-Steklov problem (see [3] for formal 
introduction of this problem). 



Ap = 

dip 
dy 

p = 



up 



in W, 
on F, 
on B. 



(2.1) 
(2.2) 
(2.3) 



Here we assume that W = W(D) is an axisymmetric liquid domain. The problem ( 2.1| 



(2.3) has a variational formulation similar to the one described in the Introduction for 



the sloshing problem (see (1.6)), with the only difference in the class of admissible 
functions for the Rayleigh quotient, which is now the space of H 1 (W) functions which 
vanish continuously at B (see [3] for more details). Since W is axisymmetric it is 
possible to use the same ansatz (1.8) for Dirichlet-Steklov problem as for the sloshing 
problem. We denote the eigenvalues of (2.1)-(2.3) by v m ,k(W) m a similar manner as for 



the sloshing problem (1.1 )— (1.4). A standard argument shows that the first eigenvalue 



of (2.1)-(2.3) is simple, it equals Vq,i(W), the corresponding eigenfunction has constant 



sign, and it is the only eigenfunction with this property. 

Using standard arguments (see e.g. [31 Section 3], [26J, [30]) one obtains the following 
domain monotonicity results for eigenvalues for both the Dirichlet-Steklov problem 



(2.1)-(2.3) and the sloshing problem (1.1 )— ( 1.4). We omit the proofs of these results 
because they are very similar to the proofs of Propositions 3.1.1 and 3.2.1 in [3]. 

Lemma 2.1. Let W\, W2 be axisymmetric liquid domains. If W\ C W2 and F± C F2 
then 5b,i(Wi) > 5b,i(W 2 ). 

Lemma 2.2. Let W\, W2 be axisymmetric liquid domains. IfW\ C W2 and F\ = F 2 
then f ,i(Wi) < ^0,1(^2) and Ui i(Wi) < ^i,i(W 2 ). 
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In the rest of this section W denotes a liquid domain satisfying assumptions of The- 
orem 



1.1 and such that F is the unit disk. By scaling it is sufficient to consider only 
such case. 

One of the important tools in the proof of Theorem (i) is the Stokes stream 



function * corresponding to ip (see [281 P- 125-127]). Suppose that ip is an axisym- 



metric eigenfunction of the sloshing problem (1.1)-( ]1.4[ ). Then * is the axisymmetric 
continuous function defined on W by 

dr r dy 1 dy r dr 



Here we use cylindrical coordinates (1.7). Since ip satisfies Neumann boundary condition 



on B, * is constant on B. Note that formula (2.4) defines * uniquely up to a constant. 
Hence we may assume that $ = on B. Furthermore, * — when r = 0, and * 
satisfies in W the relation 

<9 2 * 19* d 2 ^ 

1 = 0, r^O. 

dr 2 r dr dy 2 

In particular, * attains its maximum and minimum on {r = 0} U dW. From the 
boundary conditions it follows that the extreme values of * are attained on F . 

In this section we will need the following result which follows by a classical Courant- 
Hilbert argument (we omit the standard proof). 



Lemma 2.3. Any axisymmetric eigenfunction of the problem (1.1 )— (1.4) corresponding 
to vo,i(W) has two nodal domains. 

Lemma 2.4. Suppose that fo,i(W) < Uq i(W). Then the Stokes stream function * cor- 
responding to an axisymmetric eigenfunction ip with the eigenvalue Uo^W) has constant 
sign on W . 

Proof. Suppose, contrary to the hypothesis, that * does not have constant sign. By the 
maximum principle, * must change sign on F. Since in cylindrical coordinates (r, 6, y) 
we also have *(0, 6, 0) = *(1, 6, 0) = 0, it follows that there exist < r\ < r 2 < 1 such 
that U(n,0,O) = ™(r 2 ,6,0) = 0. But on F we have u 0>1 (W)(p =% = ~M- It follows 
that the set {ip = 0} intersects F along at least two circles r = r x and r = r 2 . Since ip 
has only two nodal domains W%, W 2 , one of them, say Wi, must touch F at the annulus 
ri < r < r 2 . Since W 2 is connected, we conclude that dW\ does not intersect B. Hence, 



tp restricted to W\ is the first eigenfunction of the spectral problem (2.1)-(2.3) in W% 
with F and B replaced by Int {y=0} (Fn<9iyi) and (dWi)\F. Hence u rfW) = i>a,i(Wi). 
Here lnt{ y=0 y(F n dW\) denotes the relative interior of the set F D dW\ in the plane 
y = 0. 

By domain monotonicity, we have Vo,i(W) = Vo,i(Wi) > v 0j i(W), a contradiction 
with the assumption Uq i(W) < Vq i(W). □ 



Lemma 2.5. Suppose that Vo t i(W) < vq,i{W). If the stream function * corresponding 

g r \ r Q r 



to an axisymmetric eigenfunction ip satisfies l-(-^)<0 on F, then ip corresponds to 



the first axisymmetric eigenvalue fo,i(VK). 
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Proof. Note that < at some point on F. Indeed, if 7r(-|r) was identically 

■> or \ r or ' 1 ' Or \ r or ' J 

zero on F, we would have * = C\ + C^r 2 on F, which contradicts * = for both r = 
and r = 1. 

Let z/ be the eigenvalue corresponding to (p. Let y?o,i be an arbitrary axisymmetric 
eigenfunction corresponding to ^o,i(WO and *o,i be its stream function. We have 



1 f dip d<p ,i 1 / 19*9* 



'0,1 



/ F ' ^^0,1 J f dy dy wo,\ Jf r2 9r dr 
Integration in polar coordinates and then integration by parts yield that 

2tt f 1 19*9*01, 2tt f 1 T 9 /19*\ , 

Wo,i = / - ^ — ^ dr = / *o,i q- —e~ dr - 

vv o,\ Jo r dr dr wq,i Jo dr \r dr J 

Since ^r(^f^) < and it is not identically zero on F, and since *o,i has constant sign 
on F (by Lemma 2.4), we obtain that f F pp ^ 7^ 0. Hence ip corresponds to z/ ,i- D 

Now we come to the key element of the proof of Theorem (i), the application 
of the inverse method from Troesch paper [35]. In [351 P- 283-284] an eigenvalue v 
corresponding to an axisymmetric eigenfunction for the sloshing problem is computed 
for a family of domains. With the notation of [33], we take 02 = 4, so that the free 
surface is a unit disc. When A G (0, 2], then for the domain 

W x = {{x,y, z):x 2 + z 2 < Ay 2 + 8y/X + 1, y < 0} 

we simply have v = A; the corresponding eigenfunction is the polynomial ip = 1 + Xy + 
4y 2 — 2r 2 + (4/3)A?/ 3 — 2Xr 2 y, where r 2 = x 2 + z 2 , and the stream function corresponding 
to p is * = (A/2)(r 2 — r 4 ) + Ar 2 y + 2Xr 2 y 2 . Since J^(-ff) < on F, we may apply 
Lemma 2.5 an we obtain the following corollary. 



Corollary 2.6. For the set W\, either ^0,1(^61) ^ ^o,i(Wa) or the Troesch eigenvalue 
v = X is equal to Vq^Wx). 

It is a natural conjecture that in fact always v = Vq i(W). However, we were not able 
to prove it. For our purposes the above corollary is enough. 

In the proof of Theorem 1.1 (i) we need some knowledge about sloshing eigenvalues 
for cylinders. These eigenvalues are well known (see e.g. [3, example 2.1, p. 24]). 
We collect some results about these eigenvalues which we need in this section in the 
following lemma. 

Lemma 2.7. Let h > 0, U h = {(x, y,z) : x 2 + z 2 < 1, —h < y < 0}, F = {(x, y, z) 



x 2 + z 2 < 1, y = 0} and B = dUh \F. Let us consider the sloshing problem (1.1) -(1.4) 
in the cylinder Uh- Then we have 

vi,i( u h) = tanh(j^/i), u ,i(U h ) = j ,i coth(j ,ih), 

where j[ 1 ~ 1.8412 is the first positive zero of J[ and j ,i ~ 2.4048 is the first positive 
zero of Jq. Here J and J\ are Bessel functions of the first kind. 
For any h> we have Vi,i{U h ) < j[ 1 and v ,i(U h ) > j 0>1 . 



ON HIGH SPOTS OF THE FUNDAMENTAL SLOSHING EIGENFUNCTIONS 



9 



Proof of Theorem 1.1 (i). Let us recall that we assume that W is a liquid domain sat- 



isfying assumptions of Theorem 1.1 and F = {(x, y, z) : x 2 + z 2 < 1, y = 0}. 



First note that any eigenfunction of (1.1)-(1.4) corresponding to ui(W) must have 2 



nodal domains. Eigenfunctions of (1.1 )— ( 1.4) corresponding to v m ,k(W), m > 2, of the 



shape il)m,k{r,y) cos(m6), ip m> k(r,y)sin(m9), have at least 4 nodal domains, so v mj k(W) 
cannot be equal to ui(W) for m > 2. Recall that vi^iW) > z/l,i(W). Hence, in order 
to show that Ui(W) = fi,i(W), we only need to prove that z/ 01 (W) > f^iiW). 

Note that for A G (0, 2] we have Ay 2 + Sy/X + l < (l + 4y/A) 2 , so that is contained 
in the circular cone 

V x = {{x, y, z):x 2 + z 2 <{l + Ay/\) 2 , y < 0}. 
For A = 2, in fact, W\ = V\. The height of V\ is equal to A/4. Hence, Ux/i is the 



smallest vertical cylinder containing V\. By Lemma 2.7 fi t i(U\/4) = Ji ]_tanh(jj iA/4) 
and Vo,x(Ux/i) = jo,i coth(j ,iA/4). 

Let h be the height of W. When h< 1/2, then for A = Ah we have W X QV\QW C 
U\u. There are two possibilities. 



(a) If vo t i(W\) > i^o,i then, by domain monotonicity and Lemma 2.7 we have 

V0,l{W X ) > V 0> i(Ux/4) > Vl,l[Ux/4). 



(b) Otherwise, by Corollary 2.6 we have ^o,i(W>) = A, and since A > (j^ 1 ) 2 A/4 > 
j^tanh^'^A^) = vi,i{Ux/a), we also have v 0>1 (Wx) > ^1,1(^/4)- 

Hence, in both cases, we have v 0s i(W) > fo,i(WA) > ^1,1(^/4) > ^1,1 (H 7 ), as desired. 

When the height h of H 7 is larger then 1/2, then we simply have W 2 = V2 C C L^- 
Again, there are two possibilities. 



(a) If ^0,1 (W2) > ^0,1(^2), then by domain monotonicity and Lemma 2.7 we have 
"o,i (W 2 ) > i>o,i (W 2 ) > u ,i{U 2 ) > j'0,1 > 

(b) Otherwise, we have ^0,1(^2) = 2 > j[ v 



As before, in both cases we have, by domain monotonicity and Lemma 2.7, vo t i(W) > 
^0,1(^2) > fit > V\,\{Uh) > Vi t i(W), which completes the proof of the inequality 

Hence Vi(W) = v\ \ (W). There are exactly 2 linearly independent eigenfunctions 
corresponding to i>i{W) = Vi t i(W): i/ii t i(r,y) cos(9), ^1,1 ( r ) y) sin(0). Hence i>i{W) has 
multiplicity 2. □ 

3. Continuous dependence under variation of the domain 

In this section we first define a new class of domains W and formulate monotonicity 
properties of ipx,i f° r this class. Then we prove continuos dependence of u^i and ipi^ 
under certain variations of the domain. Ideas used in this section are similar to the 
ideas from Section 2 in [20J. 

Let us first describe the class W in an informal way. The class W consists of axially 
symmetric domains W with horizontal cross-sections being circles or radius decreas- 
ing with depth —y. A similar condition for two-dimensional domains was assumed 
in Theorem 1.1 in |20j. For technical reasons we will assume certain regularity near 
the free surface and the vertical axis. More formally, we first describe the class D of 
cross-sections. 
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Figure 2. An example of a domain D belonging to the class D. 

Definition 3.1. A domain D C {(r, y) : r > 0, y < 0} belongs to the class of domains 
D iff its boundary consists of the following 3 parts (see Figure (2J): 

(i) the horizontal interval F(D) = {(r, y) : r G [0,ro), y = 0}, where ro > 0; 

(ii) the vertical interval R(D) = {(r,y) : r = 0, y G (jfo, 0]}, where yo < 0; 

(hi) B(D), parametrized by a simple continuous curve (r(t),y(t)), t G [0, T], satis- 
fying the following conditions: 

(a) (r(0),y(0)) = (0,y ), (r(T),y(T)) = (r ,0), and r(t) > and < are 
nondecreasing for t G [0, T], 

(b) there exist e > (e < T/2) and M > 1 such that for t e [0, e], r(t) = t and 

is a Lipschitz function with Lipschitz constant M, and for t G [T— e, T], 
y(t) = t — T and r(t) is a Lipschitz function on [T — e, T] with Lipschitz 
constant M. 

For fixed e > 0, M > 1, H > e, Tq > e, we write D G T>(e, M, H,r ) when the above 
relations hold with the prescribed e, M and r , and for some y G (— if, — e). 



Definition 3.2. The domain W C {(x,y,z) : x, z G R, y < 0} belongs to i/ie class of 
domains W iff 

W = W(D) = {(x, y, z) G E 3 : (Vx 2 + 2 2 , y) G D U -R(-D)} 
for some D E D. 

For PF = VF(-D) G W we will always assume that its boundary dW consists of 2 
nonempty parts F, B with B = dW\F and F = {(x, y, z) : x 2 + z 2 < r%, y = 0} where 
r is the number appearing in the definition of D. 

By W(e, M, H,r ) (e > 0, M > 1, H > e, r > e) we denote these domains W = 
W(D) from W for which D G D(e, M, H, r ). 



Note that all bounded, convex, axisymmetric domains satisfying the John condition 
belong to W. We are able to prove the monotonicity of ipi i stated in Theorem 1.1 for 
all W G W. 



ON HIGH SPOTS OF THE FUNDAMENTAL SLOSHING EIGENFUNCTIONS 



11 



Theorem 3.3. Let D G D, W = W{D) G W. Let be the eigenfunction of (fL9| 



(1.12) corresponding to the eigenvalue V\ \ for the domain W = W(D). After multipli 
cation by ±1 we may assume that ip\j_ > on D. Then we have 



>0 mD, (3.1) 

> in D. (3.2) 



dr 
dy 

In order to use methods from [20J we need to introduce another class of domains. 
This class may be briefly desribed as star-shaped Lipschitz domains (cf. definition of 
class L M in [12 p. 744]). 

Definition 3.4. Let p G R 3 . By <A(po) we denote the class of domains WcR 3 such 
that 

W = L G R 3 : < \p - po| < / (^^] \ U {p } 



.b — Pol 

for some positive Lipschitz function / on the unit sphere S 2 C R 3 . If < r\ < r 2 and 
rj > 0, then by A(p , r 1 ,r 2 , rj) we denote the subclass of A(p ) for which / is a Lipschitz 
function with Lipschitz constant r] such that r\ < / < r 2 . 

Lemma 3.5. Let p^ G R 3 . IfW G A(po), then W is a Lipschitz domain; that is, there 
exist constants 5>0,N£¥S,L>0 such that the boundary dW may be covered by balls 
Bi, i = 1, . . . , N of radii 8, and such that for each i, R dW is the graph of a Lipschitz 
function with Lipschitz constant L. When < r\ < r 2 , rj > and W G A(po, ri, r 2 ,r]), 
then 5, N and L depend only on n, r 2 , rj. 

The proof of this lemma is standard and is omitted. 

We will use the following notation: we will write C(a, (3, . . .) to indicate that C is a 

constant depending only on a, /3, 

The following lemma is the crucial result which is taken from [20] . 



Lemma 3.6. Let po G R 3 and let W G A(po) be a liquid domain contained in a ball 



B(p ,r 2 ). Let ip be a solution of the eigenvalue problem ( 1.1 )— ( 1.4 ) such that f F (p 2 = 1, 
and let v be the corresponding eigenvalue. Then tp G H 3 ^ 2 (W) and there is an extension 
<p G if 3//2 (R 3 ) ofcp (that is, (p = cp onW ) such thatsuppip C B(xo,r 2 +l) . In particular, 
V<p G L 3 (R 3 ). 

Suppose that < r\ < r 2 , rj > 0, u > and a > 0. If W G A(po, r%, r 2 , rj), v < u 

and \\<p\\l2(W) — a > then \\tp^H 3 / 2 (W)> ll^ll^ 3 / 2 (R 3 ) an d ||V0|U 3 (b. 3 ) are bounded above by 
constants depending only on r\, r 2 , rj, uq and a. 



Proof. The proof is based on the arguments used in the proof of Lemma 2.5 in 
For the convenience of the reader we write here much more details, but the main idea 
is exactly the same. As it will be seen below the fact that <p G H 3 ^ 2 (W) follows easily 
from [T7] and [18] . What is more difficult to justify is the fact that the norm |M|# 3 / 2 (w) 
depends only on r\, r 2 , rj, u and a, and the properties of the extension. In order to 
justify these facts we repeat some of the arguments from 
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Suppose that W G A(p , r 1; r 2 , n). Consider ip as the solution of the Neumann prob- 
lem 

dip 

Aip = on W, — = ul F ip on dW. 
on 

By M(u) we denote the nontangential maximal function of the function u : W — > R 
defined for q G dW by 

M(u)(q) = sup{|u(p)| : p G W, \p — q\ < /tdistQo, dW)} 

(cf. [TTl p. 22], or [IH])- Here we take k = 2(L + 1) where L is the constant from 
Lemma |3.5| Of course L and k depend only on r 1 , r 2 and n. By M(Vy2) we mean 

;g) + M(|) + M(g; 



Md 2 ) + M(§2) + M(ff). Theorem 4.1 in [T7] or Theorem 2 in [E] give that 



||M(V^)|U 2(aw) < CH^IIlW) = Cu 2 = Cv 2 < Cvl (3.3) 

and Corollary 5.7 in [18] (cf. also [T9l Remark (b), p. 206]) gives that ip G i/ 3/2 (H/). 

Now our aim will be to show that the norm IMI^a^^ depends only on ri, r 2 , r], u , 
a. At first we want to justify that the norm ||M(Vv2)||l2( 9H /) depends only on ri, r 2 , 
and uq. This follows from the proof of Theorem 2 in [19]. Indeed this theorem is proved 
first for star-shaped smooth domains Q = {(r, 9) : < r < ^/(6)}, \I/ G C 00 (S' n ~ 1 ), see 
[T9| p. 205, line 4]. For such domains the assertion of Theorem 2 holds with a constant 
depending only on the Lipschitz constant of the function H^Hoo; mingg^n-i |\I/(6 I )| and 
n, see [IHl p. 205, formula (5)]. The rest of the proof of Theorem 2 is the approximation 
of general star-shaped domains by smooth star-shaped domains. It occurs that the 
constant C in the formulation of Theorem 2 depends only on the Lipschitz constant of 
the function \1/ defing the star-shaped Lipschitz domain, H^Hoo, ming^s™- 1 \^(@)\ and 



k. In our situation this gives that the constant C appearing in (3.3) depends only on 
ri, r 2 and n. Note that in [19] the nontangential maximal function is defined for k = 2 
but the proof for the constant k = 2(L + 1) is exactly the same. 
Now let 

T(q) = {p eW : \p- q\ < Kdist(p,dW)}, q G dW 
and define the area integral of a function v : W — )■ R by 

1/2 

S( v )(q)=[ I dist(p, dWy^Vvip)] 2 dp' 
'r(q) 

By S(Vv) we mean 5(g) + S(^) + 5(g). By Theorem 5.4 in [18] (cf. also [TBI 



Corollary 5.7, (a) =^ (c)]) and (3.3) we obtain 

\\S(V<p)\\ L 2 {dW) < C\\M(V<f)\\ L 2 {gw) < C, 
where C = C'(ri, r 2 , n, z/ ). We claim that (cf. [HI Corollary 5.7]) 

/ dist(p, dW)\V 2 V (p)\ 2 dp < C", (3.4) 

where C" = C"(ri,r2,n, i>o). Here V 2 y? is the vector of all second derivatives of ip (cf. 
[T8| p. 181, line 6]). The inequality for the integral over a neighborhood of dW follows 
from the estimate of ^(V^)]!^^^) by appling the argument used in the proof of the 
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upper bound inequality in (6.1) of j8] to each of g, g. The integral over the 
remaining compact subset of W is bounded using a simple gradient estimate for each 
of the harmonic functions g, g, and the inequality |Vy?| 2 = z/ < u . Our claim 
is proved. 

Now we will argue like in Corollary 5.7, (c) ==>■ (b) in [18J. Note that Vy? is harmonic 
so one could use Corollary 5.5, (c) ==>• (b) in [18] for v = V<p. The implication 
(c) ==>- (b) in Corollary 5.5 in [18] follows from the proof of Theorem 4.1, (b) ==>- (a) 
in [18] for u = v = Vip. Indeed we have 

\\5 1/2 (p)Vu(p)\\ L 2 (w) + ||w|| L 2 W < C(r,R,r},is ) 

for u = Vy?, 5(p) = dist(p, dW). From the proof of Theorem 4.1, (b) ==>- (a) in [TS] it 
follows that 

u = Vy?G [L 2 (W),H 1 (W)} 1/2 , 2 

and ||V^||[l2( W ) , j h" 1 (iv)] 1/2 2 — C( r i> r 2; Vi u o)- The last inequality follows from the proof 
of Theorem 4.1, (b) =>■' (a) in [18] and Lemma |3~o] Here [L 2 (W), H\W)] W is the 



interpolation space of power 1/2 given by the real interpolation method and Vy? G 
[L 2 (W), H l (W)}i/ 2 ,2 means that g, g, g € [L 2 (H/), i/ 1 ^)] 1/2,2- 

Now we will argue in the similar way like in Proposition 2.17 in [18] (or in Proposition 
2.4 in [18] with real interpolation instead of complex interpolation). By Theorem VI. 5 
in [33] for any bounded Lipschitz domain D C R n there is an extension operator E 
mapping functions on D to functions on R n such that Ef(x) = f(x), x G D, and such 
that for each k G IN, E is a bounded operator from H k (D) to if fc (R n ) (cf. also Theorem 
2.3 in [TS]) (if = L 2 ). When _D is contained in a ball B(p,R), R > 0, the operator 

may be chosen so that for any function / on D we have supp(i?/) C B(p,R + 1). 
Assume that there exist numbers S > 0, iVG IN, L > such that the boundary dD 
may be covered by balls Bi, % = 1, . . . , N of radii 5 such that for each i, Bid dD is the 
graph of a Lipschitz function with Lipschitz constant L. It follows from Theorem VI. 5' 
in [33] and the proof of Theorem VI. 5 in [33] ( pages 190-192, see in particular (31) on 
page 191 and inequalities on page 192) that the norm of E : H h (D) —t H k (H n ) depends 
only on n, k, 8, N, L. What is more formula (1.19), page 174 from [TS], arguments 
after this formula and formula (31), page 191 in [33] give that 

— (£?/) = Q ifi (—/) + 5^/), J = l,-..,n, (3.5) 
3 i= i 3 

where Qj,i, Sj are bounded linear operators from H k (D) to H k (H n ) such that their 
norms depend only on n, k, 5, N, L. It follows that if /, V/ G H k (D) then V(Ef) G 
H k (R n ) and \\V(Ef)\\ HHMn) < C(n, k,6, N, L)(\\V f\\ HHD ) + \\f\\ HHD) ) 



Now let us come back to our situation. By Lemma 3.5 the norms of E, Qjj, Sj 
as bounded linear operators from H k (W) to H k (R 3 ) depend only on k, r 1; r 2 , 77. We 
know that tp G H\W), Vtp G [H°(W), H 1 (W)] 1 / 2 ,2 and y\\m(w) < C(n, r 2 , rj,vo, a), 



\\\/(f\\ [h° (w) ,h 1 (w)]i/ 2 2 — C(r 1 ,r 2 ,r],i i o). By the real interpolation method and (3.5) 



(cf. the proof of Proposition 2.4 in [TS]) it follows that Eip G H 1 (R 3 ), V(E<p)e 
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D' 




1 



r 



(r(a),y(a)) 



2/o + ho 
-H + ho 



Figure 3. An auxiliary picture for Lemma 3.7 



[i?°(R 3 ),i? 1 (R 3 )]i /2 , 2 and also ||^|| H i (]R 3) <C(r 1 ,r 2 ,ri,u ,a), ||V(£^)|| [H o (R 3 )iif i (]a 3 )]l/2i2 < 
C(r u r 2 ,ri,v ,a). 



in [5] and the second inclusion follows from Theorem 6.5.1 in [5]. 

Recall also that supp(Eip), supp(V(-Ey?)) C B(p ,r 2 + 1). It follows that V(Eip) G 
L 3 (B(p ,r 2 + 1)) and \\V(Eip)\\ L 3 {B{p(hr2+1)) < C(ri,r 2 ,r),is ,a). We also have Eip G 
H\B{j) Q ,r 2 + 1)), V(^) G H l ' 2 {Bip ,r 2 + 1)), and the norms H^U^^+i)), 
|| V(-E , </ , )ll_H" 1 /2(B( P0ir2+ i)) are bounded by C(ri,r 2 ,r],u ,a). It follows that ( see e.g. 
H21 formula (38)], cf. also [TH1 Proposition 2.18(a)]) E<p G H 3 / 2 {B{p ,r 2 + 1)) and 
\\E<f\\ H ^(B(p ,r 2 +i)) < C(ri,r 2 ,V,^o,a). □ 

Recall that by scaling it is sufficient to consider the domains W such that F is a unit 
disk. For this reason we only consider W G W(e, M, H, 1) (that is we fix r = 1). 

Lemma 3.7. Fix e G (0,1), M > I, H > e. Assume that D G T>(e, M, H, 1). Le£ 
fro — e/(2M) and let D' be the translation of D by the vector (0, ho) (see Figure^. Let 
(f, a) be the polar coordinate system in the (r, y) plane, with y = f cos a and r = f sin a. 
Then there is a Lip schitz junction f on [0, tt] such that D' = {0 < f < /(a), < a < 
7r}. Furthermore, the Lipschitz constant of f , the infimum of f and the supremum of 
f depend only on e, M, H . 

Proof. The proof of this lemma is standard. Let (r(t),y(t)), t G [0, T], be the curve 
from the definition of T>. Then the boundary of D' consists of a part of the ?/-axis, a 
horizontal interval, and the curve (r(t),y(t) + h ), t G [0,T]. In the polar coordinates, 
this curve is given by 



Our goal is to prove that a(t) is increasing, and that / = r oa 1 is a Lipschitz function. 



We have 




fit) = V(r(t)) 2 + W) + h y, 



a(t) — — h arctan 



y{t) + hp 
r(t) 
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By an appropriate reparametrization, with no loss of generality we assume that r(t) 
and y(t) are absolutely continuous functions of t, and r'(t) + y'(t) > for a.e. t G [0, T]. 
Since r(t) > r(e) = e for t G [e, T], and y(t) < y(T -e) = -e for t G [0, T - e], we have 
>s-h for all t. Furthermore, r'(t) = 1 and < < M for a.e. t G [0,e], and 
similarly < r'(t) < M and = 1 for a.e. t G [T - e, T]. We find that f (t) and a(t) 
are absolutely continuous, and for a.e. t G [0, T], 

^ r(t)r'(t) + + h )y'(t) r(t)y'(t) - (y(t) + h )r'(t) 

For a.e. t G [e, T — e], we have 

> gjAO + (g - Mil > e y^l+i^) 

1 J - (f(t)) 2 " 2 (f(t)) 2 ' 

For a.e. t G [0, e\, 

. = fjAg ~ (y(t) + h ) e-hp e y'{t)+r'{t) 

a[) (r(t)) 2 " (f(t)) 2 " 2(1 + M) (f(t)) 2 • 

Finally, for a.e. t G [T-£,T], 

aU " (f(t)) 2 " " - (r(t)) 2 " 2(r(t)) 2 " 2(1 + M) (f(t)) 2 ' 

In particular, a(t) is strictly increasing, with a(0) = and a(T) = n/2 + arctan(/i ). 
Furthermore, in a similar manner, we have for a.e. t G [0,T], 

. r(t)r'(t) + |y(t) + W(*) , r^+^g) 

Hence, for a.e. t G [0, T], 

I'MI < (1 + ff + e ) ?<i±M r(t) < 2 ' 1 + M " 1 + g + £ ' 2 . 



a'(t) 

It follows that £>' = {0 < f < /(a), < a < 7r} with f(a) given by f(a(t)) = r(t), 
t G [0, T] (that is, for ct G [0, &(T)]), and by /(a) = —ho/ cos a for a G [a(T),7r]. Since 
f (t) > e - h > e/2, / is bounded below by e/2. Also, /(d) < 1 + H + e. Finally, / is 
absolutely continuous, for a.e. t G [0, T] we have 



a (t) £ 



and for d > d(T), 



_ h sin d /io hl + 1 

(cosd) 2 ~~ (sin(arctan(/i ))) 2 h 
so that / is a Lipschitz function with Lipschitz constant depending on e, M, H . □ 

As an immediate conclusion of this lemma we obtain the following result. 

Corollary 3.8. Let e G (0, 1), M > I, H > e and assume that W G W(e,M,H, 1). 
Then there exist p = (0, —e/ (2M), 0) and < r\ < r 2 , rj > 1 depending only on e, M , 
H such that W G A(p , r±, r 2 , rj). 
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As in [2U] we define the distance between domains (cf. page 745 in |20j). 
Definition 3.9. Fix e G (0,1), M > 1, H > e. Put p = (0, -e/(2M), 0). Let 



Wi, W 2 G W(e, M, H, 1). Corollary [378] gives that Wi G -A(p ), so 

P 



p G R : < |p-po| < /i 



Po 



U{p } 



1,2, 



— Pol 

for some Lipschitz functions fi on the unit sphere S 2 . We define the distance between 
Wi, Wi G W(e, M, H, 1) by 

dw( e ,Af,H,l)(Wl, W 2 ) = \\fl - / 2 ||oo- 

When it is clear to which class W\, Wi belong we will abbreviate dw( e ,M,H,i)(Wi, W2) to 
d{W u W 2 ). 

Roughly speaking, we measure the distance between W±, W 2 G W(e, M, H, 1) as the 
natural distance between star-shaped domains in the class A(p ) for a special choice of 
Po depending on e, M. 

Lemma 3.10. Assume that D s G D(e, M, H, 1) for s in a neighborhood of 0, for some 
fixed e G (0, 1), M > 1, H > e. Let r s (t), y s {t) be functions defining B(D S ) in the 
sense of Definition 3.1. Assume that all functions r s , y s are defined on [0, T], and that 



r \ 



— > and 



\Vs ~ 2/0 oo 



— > as s — > 0. Then 



d(W s , W ) = d w(£M ,H,i)(W s , W ) -)• 0, as s 0, 



Proof. Let r s (t), a s (t), / s (a) be defined as in Lemma 3.7 (see Figure [3J), but for the 
domain D s . Note that f s (a) = fo(&) for a greater than a s (T) = &o(T). 

Clearly, f s (t) — > ro(t) as s — > uniformly in t G [0,T]. In a similar manner, since 
r s(£) > e for t G [e, T], we have a s (t) —¥ a (t) uniformly in t G [s,T]. Furthermore, for 
t G [0, e] we have y s (t)+h < —(e—h ) < —s/2, and a s (t) = — arctan(r s (£)/(y s (t)+/t ))- 
It follows that a s (t) — > &o(t) uniformly also in t G [0,e]. Finally, since all of f s are 
Lipschitz functions with Lipschitz constant rj, we have 

\fs(&o{t)) ~ /o(«o(t))| < \fs(a (t)) ~ fs(a s (t))\ + \f s (a s (t)) - f Q (a (t))\ 

< 7]\a (t) - a s (t)\ + \f 8 (t) -r (t)\, 

which converges to uniformly in t G [0,T]. Therefore, f s (6c) — > fo(&) uniformly in a, 
and the lemma follows. □ 



From now on we will be interested in eigenfunctions of (1.1 )— ( 1.4) corresponding to 



V\ t i. In the rest of the paper we will use the following notation: 

if>(r,y) = il>i jl (r,y), <p(r, 8, y) = V>i,i(r, y) cos 9, v = u 1A . (3.6) 

We use here cylindrical coordinates (r,9,y) defined in (1.7), cf. ( 1.8[ ). By cp(x,y,z) 
we denote the function <p(r,9,y) written in Cartesian coordinates (x,y,z). Recall 
that <f(x,y,z) is one of the eigenfunctions of ( 1.1 )—( 1.4) corresponding to v — Un, 



cp(—x,y,z) = —(p(x,y,z), and ip > on D. We will always assume that (p is normal- 
ized so that j F if 2 = 1. Since in this section we discuss the problem (1.1 )— ( 1.4 ) for more 
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than one liquid domain, we will indicate the dependence on W in the subscript, as in 

1p W , Lf W , V W . 

Lemma 3.11. Let W G W(e, M, H, 1). There exist absolute constants C\, C 2 such that 

= \Vip w \ 2 < Ci, / <£w < C 2 . 
Jw Jw 

Proof. Note that W is a subset of a cylinder {(x, y, z) : x 2 + z 2 < 1, — iJ < y < 0}. By 
Lemmas 



2.2 



and 2.7 we get ^ < j( x . By the variational characterization of vy/ and 

I F = 1 we have *V = Jvf I V V?w| 2 - 

Let H4 be the orthogonal projection of W on the yz-pl&ne. For any p G H4, let 
be the cross-section of W with the line parallel to the x-axis and passing through p. 
Then l(p) is an interval symmetric with respect to the yz-plane, and l(p) has length 
\l(p)\ < 2. Since ^wO^j Vi z ) is °dd with respect to the x-variable, for any p G W x we 
have 

^ < ML 2 / ( 9( fw\ 2 < ± r (&pw 

hp) W ~ 7T 2 J up) V 9x J ~ n 2 J Kp) V dx 



Hence 




( ftyw X 
\ dx J 



dp 



( d<p w \ 



2 4 



The following lemma is analogous to Lemma 2.5 in 

Lemma 3.12. There is a constant C{e,M,H) such that if W l ,W 2 G W(e,M,H,l) 
and vwi, v w 2 are corresponding eigenvalues v for domains W\, W 2 then 

\vw x ~vw % \< C(e,M,H)d(W 1 ,W 2 ) 1/3 . 



Proof. By Corollary |3j[ Wt,W 2 G A(p , r h r 2 , rf) for p = (0, -e/(2M), 0) G R 3 and 
some < r i < r 2 and 77 depending only on e, M, H. Using this, Lemma |3.11| and 
Lemma 3.6 we obtain that <fi G H 3 ' 2 (Wi) and ll^w^llffs/sfwi) < C\ — Ci(e,M,H) for 
i — 1, 2! Also by Lemma |3.6 there exist extensions <pyy. (i = 1, 2) of functions 
such that supp(^vvj C B(p ,r 2 + 1) and [|^w<|| ^3/2^35 < C 2 = C 2 (e,M,H). By a 
symmetrization argument, we may also assume that these extensions (pwi are again 
odd functions of x. Lemma 3.6 gives also that || V(^m/J|l 3 (r 3 ) < C3 = C^e, M, H). 
We have 

/ |V^ 2 | 2 = f \Vif W2 \ 2 + [ \V(pw 2 \ 2 -l |V^ 2 | 2 . (3.7) 
JWi Jw 2 Jw 1 \w 2 Jw 2 \W! 

By Holder's inequality (for p = 3, q = 3/2) we get 

/ |Vy5w 2 | 2= / 1vki\w 2 I Ww 2 | 2 < ||lwi\w 2 IU 3 (H 3 )ll |V^w 2 | 2 |L 3 / 2 (r 3 ) 
= \W X \ W 2 \ l/3 \\ V<pw 2 lli 3 (R 3 ) < c 4 \w x \ W 2 \ 1/3 < C b d{W 1 , W 2 ) l! \ 
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where C 4 , C 5 are constants depending only on e, M and H. 
Note also that f Wa \V<fw^ = »w 2 and j F (p 2 w% = j F <fw 2 



1. By the variational 



characterization of v and (3.7) it follows that 



fw, \ V( Pw 2 



If^w 2 ( x ^^ z ) dx dz 



<isw 2 +C 5 d(W u W 2 ) 1/3 . 



In the similar way we get i>w 2 < vw\ + Csd(Wi, W2) 1 ■ 



□ 



Let W G W(e, M, H, 1). Recall that the eigenfunction ip w has the lowest eigenvalue 
i>w among all the eigenfunctions of the problem (1.1 - 1.4[ ) on W which are odd functions 
of x, and the next such eigenvalue is strictly greater. We denote it by v w + e w , where 
Ew > is a constant depending on W. Hence if g G ^(W), j F ipwg = and g is odd 
with respect to x, then 



W JF 



(3.8) 



Lemma 3.13. Let W u W 2 G W(e, M, H, 1). There is a constant C = C(e, e Wl ,M, H) 
such that 



J \<p Wl - ip W2 \ 2 < Cd{W 1) W 2 ) l '\ 
1 > J mm > 1 - ^d(W 1 , W 2 ) l '\ 



(3.9) 
(3.10) 



Proof. Recall that ipw, {i = 1, 2) are normalized so that J F Iv^wJ 2 — 1> fWi > on W + 
and ip W - < on W-, where W + = {(x,y,z) e W : x > 0}, W_ = z) e W : x < 

0}. 

First note that (3.10) follows easily from ( |3.9 ) and the equality 



(m - m) 2 + / v>w 1 + vw 2 

J F J F 



Hence, it is sufficient to show (3.9). Put a\ = fp(pwi<Pw 2 - By the normalization of 
ipw! and <fw 2 we have < a% < 1. Let g = <£>vk 2 — ^i^Wd where <£w 2 are the 
extensions of <^wn (pw 2 defined in the proof of Lemma 3.12 On W\ D W 2 , (pw l — ¥w t 
are continuous (i = 1,2), and so also g is continuous on W\ D W^. In particular g is 
continuous on F. We also have g G L 2 (F), and by the definition of «i we have 



<p Wl 9 = 0. 



(3.11) 
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Furthermore, g G H 3 I 2 {W\) C H 1 (Wi). By arguments as in the proof of Lemma 



3.12 



h>W 2 



w 2 \ 



u 2 



|V^ 2 | 2 + / \V<fW 2 \ 
Wi JW 2 \Wi 



|Vy5 



w 2 \ 



Wi\W 2 



(3.12) 



> / \^'Pw 2 \ 2 -C 1 d{Wi,W 2 ) 1 ^ 
= [ \V(a lVWl +g)\ 2 -C 1 d(W 1 ,W 2 ) 1 / 3 

JWi 

= a\! \V<p Wl \ 2 + 2a l [ V^ Wl Vg+ f \Vg\ 2 - C 1 d(W 1 ,W 2 ) 1 / 3 , 

JWi JWi JWi 

where C\ = Ci(e, M, H). We have f w \VtpwA 2 = vw x -> anc ^ by the Green's formula and 

'dip Wl 



(3.11 ) we get 



Hi 



= vw x \ <fwi9 = 0. 
Jf 

Since g is odd with respect to x (because ipw 1 and (pw 2 are odd) and g G if 1 (W / i), by 
d3lB and (ESb we get 



|V#| 2 > {v Wl +e Wl ) / g 



u i 



Put «2 = (f F 9 2 y^ 2 - By the definition of «i and the normalization of tpw 2 we have 
a 2 + a\ = 1. From (3.12) we get 

i/ W2 > + a^i^ + ewi) - Cid(Wi, W 2 ) l/3 

= v Wl +a 2 2 e Wl -C 1 d(W 1 ,W 2 ) 1 / 3 , 



where C\ = Cx(e, M, H). Using this and Lemma 3.12 we obtain 



ai < 



m 2 - v Wl + c 1 d(w 1 , w 2 y/ 3 ^ c 2 d{w 1 , w 2 y 



/3 



< 



ew 1 



where C 2 = C 2 (e, M, H). Finally we have 

lm - <fiw 2 \ 2 = 2 - 2ax < 2 - 2a{ = 2a z 2 < 
Here we used the fact that < ol\ < 1. 



2 o o„ o„2 o 2^2C7 2 rf(H/ 1 ,H/ 2 ) 1 /3 



□ 



Lemma 3.14. Let W u W 2 G W(e, M, H, 1). T/jere zs a constant C = C(e, e Wl ,M, H) 
such that 



I \V( VWl -^ w . 2 )\ 2 < Cd(Wi,W 2 )^ 3 . 

JWiDWo 



(3.13) 



Proof. We have 



Winw 2 



111 



111 



? w 2 | - 2 / Vy^Vy^. 

'Wi 



3.12 
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Note that f w \V<pwi\ 2 — v w x an d, as in the proof of Lemma 

|V<^y 2 | 2 = / |V^ 2 | 2 -/ |V^ 2 | 2 + / |V0w 2 | 
Wi JW 2 JW 2 \Wi JWi\W 2 

< u W2 +C 1 d(W 1 ,W 2 ) 1/3 

< v Wl +C 2 d(W 1 ,W 2 ) 1/3 



where C\ = C±(e, M, H) and C 2 = C 2 {e 1 M, H). By the Green's formula and Lemma 3.13 
we obtain 

j V Wl V^, = j W , + j F (^) W , - Jj^wMw, 



JF 

>v Wl -^d(W u W 2 )Vz 

>p W L-C 4 d(W h W 2 ) 1/3 , 

where C?„ = C^e.ewn M, H) and C4 = C±(e, Ew 1 , M, H). The last inequality follows 
from Lemma [3.11| Finally, we obtain 



/ |V(pwi - <Pw 2 )\ 2 < 2u Wl + C 2 d(W 1} W 2 ) 1/3 - 2u Wl + 2C 4 d(W u W 2 f 



1/3 

I v \fWi fW 2 )\ ^ ""W\ I ^A^\ ' ' Li " A) ^"W\ I "^4"\ ' ' Li ' ' A} 

, w 1 nw 2 

which gives the assertion of the lemma. □ 



Lemma 3.15. Let W u W 2 G W(e, M, H, 1). There is a constant C = C(e, e Wl , M, H) 
such that 

! \ ( p Wl - (fW2 \ 2 <Cd(W 1 ,W 2 ) 1 / 3 . (3.14) 

Proof. Note that VFiflW^ is symmetric with respect to the yz-plane, and its intersections 
with lines parallel to the x-axis are intervals. Since cpwn { Pw 2 are odd with respect to 
the x-variable, as in the proof of Lemma 3.11 we get 

/ \¥>Wl ~ ¥w 2 \ 2 < \ [ \V{v>Wi ~¥w 2 )\ 2 - 
JWinw 2 k JWinw 2 

Now the assertion of the lemma follows from Lemma 13.141 □ 

The next lemma is analogous to Lemma 2.6 in |20j . 

Lemma 3.16. Let Wt € W(e,M,H, 1) for t in a neighborhood of 0, and suppose that 
d{Wti Wo) — > as t — > 0. Let K be a compact subset of Wo- Then there exists 5 > 
such that for all \t\ < 5 we have 

K C W t and \\(<Pw t ~ <Pw )1k\\oo — > as t — > 0. 

Proof. Since K is compact we have dist(if, dW ) = 2r for some r > 0. Since d(W t , Wo) 
there exists 5 > such that for all \t\ < 5 we have K C W t and dist(-fT, dW t ) > r. 
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Figure 4. (a) An example of a domain D in the class Di. (b) An 
example of a domain W in the class W x . 

Let \t\ < 5 and p G K. Then B(p,r) C W t H W and <fw t i fw are harmonic in 
B(p,r). Using this and Lemma 3.15 we get 

WwAp) ~ <Pw (p)\ < i B / v. / \<Pw t ~ fw \ 

\B(p,r)\ J B 



< 



1 

V\ B (P> r )\ 
1 



1/2 



1/2 



'B(p,r) 



w t nw 



< C(e,e Woi M,H)r-^ 2 d(W u Wo) 1/6 , 
which implies the assertion of the lemma. □ 

4. MONOTONICITY OF THE ODD EIGENFUNCTION FOR SOME CLASS OF PIECEWISE 

SMOOTH DOMAINS 



This section is similar to Section 3 in [20J. The aim of this section is to show The- 
orem |3.3| for domains in some special subclass Wi of the class W, defined below (see 
Figure |4j(b)). First we need to define the subclass T>i of the class D. In the whole 
section we use notation (3.6). 

Definition 4.1. Let y < 0, r x e (0, 1), T = n - y , n = T x < T 2 < T 3 < T. The 
domain D C {(r, y) : r > 0, y < 0} belongs to the class of domains T)\ iff its boundary 
consists of the following 3 parts (see Figure Hj^a)): 

(i) the horizontal interval F(D) = {(r, y) : r € [0, 1), y = 0}; 

(ii) the vertical interval R(D) = {(r, y) : r = 0, y G (yo,0}}; 

(iii) B(D), parametrized by a simple continuous curve (r(t),y(t)), t G [0, T], satis- 
fying the following conditions: 

(a) (r(t),y(t)) = (t,y ) for t e [0,^], 

(b) (r(t),y(t)) = (r^yo + t-Tt) forte [T 1; T 2 ], 
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(c) y{f) = yo + t — Ti for t G [T 2 ,T 3 ], r(t) is a strictly increasing function on 

[T2 T3] 

(d) (r(t),y\t)) = (l,i/o + *- 7i) for t G [T 3 ,T], 

(e) r(t) is a C°° function on [Ti,T]. 

We denote by B X (D), B 2 (D), B 3 (D), B 4 (D), the parts of B(D) corresponding to t G 
[T 3 ,T), t G (T 2 ,T 3 ), t G (T 1; T 2 ], t G (0,T\] respectively. 

Definition 4.2. The domain W belongs to the class of domains Wi iff W = W(D) for 
some D G T) 1 . For W = W{D) G W 1 we denote 5; = {(x,y,z) G R 3 : (V2: 2 + z 2 ,y) G 
Bi(D)}, i = 1,2, 3,4. 

Note that for £ = min(ri,T — T3), 77 = — yo + 1, the domain D G Di belongs to the 
class D(e, 1, if, 1) and G W(e, 1, if, 1). 

First we need the following general lemma. 

Lemma 4.3. Let R > 0, (x, y) be a rectangular coordinate system in R 2 and -6(0, R) = 
{{x,y) : x 2 + y 2 < R 2 }. Let f G C 2 '\-R,R) be such that f(0) = f'(0) = and 
r(/) = {(x,f(x)) : x G (— R, R)} H B(0,R) be the part of the graph of f contained in 
B(0,R). Assume that g e C 2 ' l (B(0,R)) and V#(0,0) = 0. 

V IrO^y) = / or a ^ (^jZ/) e r(/), where -j^ is the normal derivative to T(f) at 
point (x,y) G r(/) ; then 

d 2 q , / \ 

Proof. By assumptions on / there exists r > such that if x G (— r, r) then (x, f(x)) G 
-8(0, i?) so (x, f{x)) G r(/). In the whole proof we will assume that x G (— r, r). 
The unit (upper) normal derivative to r(/) at (x, f(x)) G T(/) is equal to 

-f Wttt + 



dn y/l + (f (x)) 2 V <9x <9y 

The condition ^-(x,y) = gives 

|?(x,/(x)) = /'(x)||(x,/(x)). (4.2) 

By our assumptions on / we get f(x) = 0(x 2 ) and f'(x) = 0(\x\). Simil arly, as- 
sumptions on g give %{x,y) = 0(a/x 2 + y 2 ) so ^(x,f(x)) = 0(\x\). By (42) we 
get 

^(x,/(a;)) = 0(a; 2 ). (4.3) 



By Taylor expansion for ^| we get for (x, y) G 5(0, i?) 



It follows that 



= £§ (0 ' 0)f + W 2 ^ 0)m + o(|£|) - 
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By ( pL3| ) and /(£) = 0(x 2 ) we get 

d 2 g 



dxdy 



(0,0)x = o(\x\). 



This implies (4.1). □ 

Lemma 4.4. Lei D £ r D 1 and W = W(D) £ W\. We have p> £ C^iW). 

Proof. This is a standard result. The lemma follows from [221 formula (13.14)] and [2UI 
p. 63, lines 1-3], see also [32]. D 

As an immediate conlusion of this lemma we obtain: 

Corollary 4.5. Let D £ V 1 and W = W{D) £ W x . Then §£, are bounded on W. 

Recall our notation W + = {(x, y,z) £ W : x > 0}. Recall also that we may assume 
that ip > on D. 

Lemma 4.6. Let W = W(D) be an axisymmetric liquid domain. If |^ > onD and 
f£ > on D then ^ > on W+ and ^ > on W . 

ay — ay — T ox — 



Proof. In the cylindrical coordinates (r, 9, y) (see (1.7)) we have W+ = {9 £ (— tt/2, tt/2)}. 
Clearly, 

dp dip 

- = -(r,y) cos 9. 
Hence ^ > on D implies |^ > on W+. Furthermore, 



dp> dp dp ( — sin 9 

dx dr d9 



Since p = ip(r,y) cos 9, we obtain 



*-*fc|).rf# + *r..)f— )• (4.4) 



dx dr 

It follows that ip > and f£ > on D implies |e > on W. □ 



Note that if D G Di and W = W(Z?) e Wi then, in view of Lemma 4.4 , ( |44| ) holds 
on W 7 , given that r = \/ x 2 + z 2 ^ 0. 

Lemma 4.7. Let D £ D 1 and W = W(D) £ Wi. Assume that (r,y) £ D. For x = r 
and any n, m £ IN we have 

dx n dy m ^ ,V,0>) = dr n dy^ r,V ^ 
Proof. This follows from the fact that for x = r we have p(x,y,0) = ip(r, y). □ 

The proof of the following lemma uses some ideas from the proof of Proposition 3.2 
in M\. 



Lemma 4.8. Let D £ Dj and W = W{D) £ Wi. If ^ > on D and ^ > on D 
then §J > on Bi(D) U B 2 (D) U B 3 (D) (see Figure^a)). 
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Proof. Since g > on D we have > on 51(D) U B 2 (D) U 53(D) so we only need 
to exclude the possibility that ^ = at some point of Bi(D) U B 2 {D) U B 3 (D). 

On the contrary assume that §^(r*, y*) = at some point (r*,y#) G B^D) UB 2 (D) U 
B 3 (D). Since y*) = 0, -^{r*,y*) = and since the normal vector is not parallel to 

the y-axis, we have Vip(n,y*) = 0. Let T = (a, /3) G R 2 , a > 0, (5 > 0, a 2 + (5 2 = 1 be 
the tangent unit vector to B(D) at point (r*, ?/*) and iV = (— /3, a) G R 2 be the normal 
unit vector to B(D) at point (r*,y*) (see Figure 5). Let us introduce a rectangular 
coordinate system (x,y) with origin at (r*,y*) so that x-axis is in the direction T 
and y-axis is in the direction N. Note that in these new coordinates B(D) in some 
ball centered in the origin is the graph of a C 2,1 (even C°°) function / such that 
/(0) = /'(0) = 0. 



By Lemma 4.4 tp G C 2 ' l (D). Using this by [TU Lemma 6.37] there is a C 2,1 extension 
of the function ip to some ball centered at (r*, y*). We will denote this extension by the 
same letter ip. 



Then the assumptions of Lemma 4.3 are satisfied which gives 



d 2 %b 

-(tv,j/*) = 0. (4.5) 



dxdy 



Since V^(r*,y*) = we get §f (r*, y*) = 0. 

Note that = a^- + and define w = + fij?-. By Lemma 4.7, u(x,y,0) 



dx dr dy dx ' dy ' 

§§(?", y) (with r — x). Note that u is harmonic in W and by Lemma 
that a > 0, (3 > we have it > in W+. Finally, for x* = r*, 

rt(x*,y*,0) = "Tpr (?"* 5 y* ) = 
ox 



4.6 



and the fact 
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and 

du d 2 ib 

At the point (x*,y*,0) G dW+ the domain W + satisfies the inner ball condition. Hence 
(4.6) gives contradiction with the Hopf's lemma for the harmonic function u. □ 

Lemma 4.9. Let u be a function which is harmonic in a bounded domain Q C R n 
and continuous in Q. Let Qi C dVt be such that for each p G Q\ the boundary dVt 
has a tangent plane at p, the outer normal derivative §^(p) exists and §^(p) = 0. Let 
Q 2 = dQ \Q\. If u > on Q2 then u > on Q. 

Proof. On the contrary assume that there exists po G Q such that u(po) < 0. By the 
maximum principle for u there exists p\ G Q\ such that > u{p\) = min g ^w(p). 
By the normal derivative lemma (see e.g. [12] Lemma 2.33) §^(pi) < which gives a 
contradiction. □ 

Let us denote R^ = {(x,y,z) : x > 0}, B+ = BnR 3 + , F+ = FnR 3 + , B i+ = £?i DR+, 



1,2,3,4. 

Lemma 4.10. Let D e r D l and W = W(D) eW x . If g > on B 2 (D) then g > 



on D and^>0 on W+. 

dy — + 

Proof. We have 



_ = _(r,y)ooB6l 
M->n™ n„( n\ ™ ro ^ 



on W, so by the assumption ^ > on B 2 (D) we get ^£ > on S 2 +- 

Let us denote u = Then w is harmonic in and continuous in W+. We have 

m > on B 2+ . We also have « = on B i+ . Since </? = on {(x, y,z) E W : x = 0} we 
get m = on {(x, 2) G : x = 0}. Recall that ip > on D so > on W+. We get 
m = |2 = z/<^ > on F+, so m > on ~F\. Note also that ^ = on 5i + U .63+. Now the 
assertion of the lemma follows from Lemma 4.9 for Q = W + and Q\ = Bi + U -B3+. □ 

Lemma 4.11. Let D G D(e,M,H,l) and W = W(D) G W(e,M,H,l). Assume that 

1 /2 

u G m zs odd with respect to the x-axis, a = (J F ii 2 ) > and | J p (pwu\ 7^ 

Then we have 



a. 



|Vw| 2 > v w I u 2 . 



IV 



Proof. Put (3 = f p ip w u and g = w— /3y?iy. We have J F <^wg = and J F g 2 = a 2 —(3 2 ^ 0. 
We also have 

/ \Vu\ 2 = [ \Vg\ 2 + [ \(3Vip w \ 2 + 2 [ VgV^ w . 
Jw Jw Jw Jw 



By Green's formula one easily get j w VgV(fw — (cf. proof of Lemma 3.13). We have 
Jw l^fwl 2 — v w- Using (3.8) we obtain f w \Vg\ 2 > v w (a 2 — fi 2 ). This implies the 
assertion of the lemma. □ 
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Figure 6. The set V+. 



Lemma 4.12. Let D eV 1 and W = W(D) G W a 
on W. 



// g > on B 2 (D) then 



¥ > 

ox — 



Proof. Since B 2 (D) is a graph of an increasing function y 



9(r), & 



g > on B 2 (D), we get §£ > on B 2 (D). We also have 



so g > on Bi(D) U B 2 (D) U B 3 (D). By (|4.4|) 



W_ 



9r 



. ; - on B 2 (D) and 
= on Bi(D) U 53(D) 

> 0on^Ui^U4 Let 
{(x, eW: x < 0} and W = {(x,y, J7~G VF : x = 0}. Since < on W_ 



dx 



and ip > on W + we get ^ > on W . 

Assume, contrary to the hypothesis of the lemma, that there exists (x*, y*, z*) G W 
such that §^(x=k, ?/*, £*) < 0. Since ^(x,y,z) = ^(—x,y,z) we may assume that 
) G W + . Let V+ be the connected component of the set {(x, y, z) G W : 
^(x,y,z) < 0} containing (x*,y*,z*). Since || > on ^! U 5 2 U 5 3 U VFo we have 
^+ C 

Let 

F(V+) = lnt {y=0 }(9V + n F) and -B 4 (V+) = Int {2/=2/o} (dy + H 5 4 ) 

(see Figure pi). By Int{ y=0 }, Int{ y=yo }, we denote the relative interior of a set in a 
2-dimensional plane. Recall that —yo is the depth of a liquid domain D G T>i, see 
Definition |4~T1 

Since |^ > on £>! U £> 2 U B 3 U Wo and is the connected component of the set 
{(x, y,z)eW: g(x, y, z) < 0} we have g = on <9y + \ U S 4 (F+)). We also 

have 

9 / dtp" 
dn V <9x 







Oil 



B±{V+ 



and 



<9n 



dtp 
dx 



d_ 
dy 



dip 
dx 



v w 



dip 
dx 



Oil 



Let V_ = {(x,y,2;) G PF : (—x,y,z) G Of course VI C VF_. Let us define the 

function u in the following way. We put u — on V+, u — — jjj^ on V- and u = on 
VF \ (V+ U VI). We have u G C(W) and w is odd with respect to the x-axis. By Green's 
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formula we also get 



|Vd 2 = 2 / |Vw| 2 



w Jv+ 

du f du f du f . . . \ 

JT U+ / JT U+ / 7T U - ( A «H ( 4 - 7 ) 

F(V+) an JB4V+) ° n JdV + \(F(V + )UB 4 (V + )) an JV+ J 

= 2v w / u 2 = u w / u 2 . 

Jf(v+) Jf 

We have F(V + ) ^ and j F u 2 > because otherwise u = on V+ which is impossible. 
Of course, there are points on lnt{ y=Q }(dW + PI F) for which || > 0, so i^(V+) 7^ 
Int{j, = o} fl F). This implies that tt and ipw are not linearly dependent on F, so 



j p Lp w u\ ^ I f^w 2 ) 1 / 2 . By Lemma 4.11 we get 



|Vw| 2 > vw I u 2 , 



which contradicts (4.7). □ 



Lemma 4.13. Let D e Di and W = W(D) G Wi. Assume that there exists h > 
such that W contains the cylinder {(x,y,z) : x 2 + z 2 < 1, > y > —h}. Then there 
exists C(h) such that for all (x,y,z) G W such that y < —h we have 

\<p(x,y,z)\<C(h). 

Proof. Denote by W the set {(x,y,z) G W : —h/2 > y}. Let M be the supremum of 
\<p\ over W . By the maximum principle, the supremum is attained at the boundary of 
W . By the normal derivative lemma (see [121 Lemma 2.33]), since ip is not constant, it 
cannot attain its supremum or infimum at dW fl {(x, y, z) : y < —h}, the part of the 
boundary where ip satisfies Neumann boundary condition. It follows that M = \ip(p )\ 
for some p = (x , —h/2, z ) G dW . 

By the Harnack inequality up to the Neumann boundary (see [U Theorem 3.9]), there 
is 5 = 5(h) G (0, h/2 A 1) such that <p(p) > M/2 for p G B(p 0l 5) fl W. Furthermore, 
\B(po,5) fl W\ > Ci5 3 , where C\ is an absolute constant. It follows that 

, f M 2 d5 3 M 2 

<P > > — • 

w JB(p th 5)nw 4 4 



On the other hand, by Lemma 3.11 , the left-hand side is bounded above by an absolute 



constant. □ 

For a fixed domain D G Di we define the family of domains D s , s G [0,1] such that 
Di = D and D is a rectangle (0, 1) x (y , 0) (see Figure [7]). Let y , r 1; r(t), y(t) be as 
in the definition of class T)\ for the domain D. For y G [yo, 0] put g(y) = r{r\ +yo — y). 
Note that D = {(r,y) : y G {y ,0), r G (0,g(y))}. The domain D s G t D 1 for s G (0, 1) 
is defined by 

D s = {{r,y) : y G (y ,0), r G (0, 1 - s + sg{y))}. 
We denote 1 - s + sg(y) = g 8 (y), so that D s = {(r,y) : y G (y Q ,0), r G (0,g a (y))}- 



2<s 
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Figure 7. Family of domains D s 
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D J 










R 


Q 



Figure 8. An auxiliary picture for Lemma 4.14 



Lemma 4.14. Let D G D\. We define a rectangle R (see Figure^ by 

i2=(n/2 J 3/2)x( 3 /f,^) J 

wherey^ = yo + fa-Tx)^, y% = -(T-T 3 )/2 andr x , y , T\, T 2 , T 3 , T are taken from 
the definition of class T)\ for the domain D. Let D s , s G [0, 1], be the family of domains 
constructed before the statement of this lemma for the domain D, and let W s = W(D S ). 
There exist a constant a = a(D) such that for all s G [0, 1] 

Qa+P 



sup 



E 



Qa r Q/3 y 



< a. 



P eD s nR a+f}=2 

Moreover, there exist constants e = e(D), H = H(D) such that for all s, q G [0, 1] 

W s G W(e, 1, H, 1) and d(W s , W q ) -> as q. 

The key point in this lemma is that the constant a does not depend on s. The proof 
of this result follows by elementary but tedious calculations. 



Proof. We define a rectangle (see Figure [8| Q = (ri/4, 2) x (yf^f), where y® = 
Ho + (T 2 — Ti)/A and y® = — (T — T 3 )/4. For each s G [0, 1] let us consider the following 
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change of variables: 

r 

y = y, r = 

9s{y) 

This change of variables is chosen so that the curve Q PI dD s is straightened in these 
new variables. In fact this transformation of variables changes Q PI dD s to {(f,y): y G 
(y?,y?), f = 1}. 

Since # s (y) = ^(y) = sg'{y), and ^(^y) = ~gf^ = we have 

did d d rsg'{y) d d 2 Id 2 



dr g s {y) dr' dy dy 9s{y) dr' dr 2 (g s (y)) 2 dr 2 ' 

d 2 d 2 2rsg'(y) cP_ + / fsg'(y) \ 2 d 2 + 2fs 2 (#'(i/)) 2 - fsg"{y)g s {y) d 



dy 2 dy 2 g s (y) dydr \ g s (y) J dr 2 (9s(y)) 2 dr' 

The set Q fl D s is transformed in the new variables to the set Q s := Ur,y) : y G 



j 2/2)3 r e ( r i/(4(7 s (y)), 1)}. Let us put ^ = t/^. Since ip s satisfies (1.9) in D s for 
m = 1, we obtain that, in the new coordinates, ip s satisfies in Q s : 

1 + {fsg'(y)) 2 d 2 ip s d 2 ip s 2rsg'{y) d 2 ip s 



(9s(y)) 2 dr 2 dy 2 g s (y) dydf 

| l + 2(rsg'(y)) 2 -f 2 sg"(y)g s (y) dij s 1 Q [ ' } 

f(g s (y)) 2 dr (fg s (y)) 2 

Now we will verify the assumptions of [T4"l Lemma 6.29]. Note that < r\ < g s {y) < 1, 
and < rsg'(y)/g s (y) < C\(D) on Q s . Hence, by an elementary calculation, there 
exists a constant A = \(D) > such that for all £, 77 G R, s G [0, 1], (f, y) G we have 
(f ~ {rsg'{y)/g s {y))r]) 2 + {l/{g s {y)) 2 W > K^ + V 2 )- This means that for all s G [0, 1], 



the operator in (4.9) is strictly elliptic on fi s with a constant A which does not depend 
on s. 

Note that the unit outer normal derivative S- on Q fl D s equals 

d 1 (d d_ 

dn {l + {sg'{y)) 2 y/ 2 \dr S9{V) dy, 

Since ip s satisfies ^ = on Q fl dD s we obtain that in the new coordinates ip s satisfies 

l + {sg'{y)) 2 fdip s l( ^dil) s 

g s (y) ~df {r ' V) ~ 39 iy) W (r ' y) = °' ( 0) 

for f = 1 and y G (yf^f)- Note that the coefficient at ^f- in the above formula is 
bounded from below by k — 1, which does not depend on s. 

We will use [TU Lemma 6.29] for a — 1, u — ip s , Q = Q s , T — {(r,y) : r — 1, y G 
(2/1,2/2)} an d ^ the operator in (4.9). We have already checked that k > does not 



depend on s and that the operator L on fi s is strictly elliptic with a constant A not 
depending on s. It is also easy to check that absolute values of all the coefficients of 



L on f2 s and the coefficients in (4.10) on T are bounded from above by a constant not 
depending on s. 
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Note that R fl D s is transformed in new variables to the set U s := {(f, y) : y G 
(yf,yf), f G (ri/(2^ a (j|)), 1)}. It is easy to check that 



sup 



d a rdPy 



ips{r,y) 



< C{D) sup_ 

(f,y)EU a a+ p< 2 



d a rd$y 



-i>s{r,y) 



(r,y)<=D s nR a+/3=2 

Note that dist({7 s , dfl s \T) is bounded from below by a positive constant not depending 
on s. By [T4"l Lemma 6.29] for u, Q, T and L as above we get 



sup_ 



d a fdPy 



<C(D) sup_|^(r, 



We have sup ( ^ )eU7 \ip s {r,y)\ = sup (rjJ/)gS7S B Y Lemma 

by a constant not depending on s. This implies (4.8). 



4.13 



this is bounded 



It is clear that for e = min(n/2, T - T 3 ), H = -y + 1 we have W a G W (e, 1, fl", 1) 

□ 



for any s G [0, 1]. The convergence d(W s , W q ) — > follows by Lemma 3.10 



en 



Lemma 4.15. Let D G T>\. Let D s , W s , s G [0, 1], be as in Lemma 4-^4\ ( see Figure 
[?|). Fix q G [0, 1]. Let R be the rectangle defined in Lemma 4-14 ( see FigureUfy. Th 
we have 

1 d 



sup 



p£D s nD q nR 



—7-0 as s — > q. 



for all s G [0, 1], and d(W s , W q ) — > as s — > q. Hence, by Lemma 3.14 

dip s dip q 2 



Proof. Denote ip s = ipD s - By Lemma 4.14, there are e, H such that W s G W(e, 1,H,1) 
a,n< 



'D,nD n 



dy dy 



->■ 



as s — >■ g. 



(4.11) 



Again by Lemma 



4.14 



^ is a Lipschitz function on D s fl i? with Lipschitz constant 
a = a(Z)) not depending on s G [0, 1]. Let 5 G (0, 1). It follows that for any p,p' G 
D s H D q C\ R with \p — p'\ < 5 we have 



9y 



(p) 



c 

9V. 



(P) 



< 



d^ t/JS dip q , 



\p' 



(p 



dy dy 

Furthermore, \B(p, 6) H D S H D q (l R\ > C(D)5 2 . Hence, 

2 



2ac^ 1 <S^' J + 2 

2 



9y 



(p 



(p) 



d Sv) 



< 8a 2 5 2 + 



<9|/ dy 
Since r > rj/2 on R, we conclude that 

2 



C(D)<5 2 



/ 




dip q 


J B( P ,6)nD a nD q nR 


dy 


dy 



sup 

peD a nD q nR 



dip s ,_ s dip q 



dy 



(P) 



dy 



(P) 



< 8a 2 6 2 + 



C(D)6 2 n 



D s nD„ 



dtps dip q 



dy dy 



By ( |4.11[ ), it follows that 
lim sup 



sup 

p€D s nD q nR 



dip s/ ^ di/j q 



dy 



(P) 



(P) 



< 2^05. 



Since 5 G (0, 1) was arbitrary, the proof is complete. 



□ 



ON HIGH SPOTS OF THE FUNDAMENTAL SLOSHING EIGENFUNCTIONS 



31 



Lemma 4.16. Let h > and W = {(x, y,z) : x 2 + z 2 < 1, > y > —h} be a cylinder. 



Then ^ > on W+. In particular 



g>0 onD = {(x,y): x E (0, 1), y E {-h, 0)}. 



Proof. The result follows from explicit formulas for sloshing eigenf unctions in cylinders, 
we have ipw( r ,y) = c (h)J\(j[ \r) cosh(jJ + where Ji is the Bessel function of the 
first kind, and j[ 1 is the first zero of its derivative (see e.g. [261 P a g e 502] or [3> page 
24]). □ 



The next result shows that the assertion of Theorem |3.3| holds for W E The 
proof of this result uses some ideas from Proposition 3.2 [20] . 

Proposition 4.17. Let D eV 1 and W = W(D) E W l . Then we have > on W 
and^>0onW+. 

Proof. In view of Lemmas 



4.10 



and 



4.12 



it is enough to show that > on D. On 



the contrary assume that there exists a domain D E T>i and a point p E D such that 

g(p)<o. 

Let D s , s E [0, 1] be the family of domains constructed before the statement of 
Lemma 4.14 for the above domain D = D\, and denote ip s = ipD a , <Ps 



Let 



q = inf I s E [0, 1] : ~w—(p) < for some p E D t 



We will first show that 



dy 



9% >Q 



dy 



on D„ 



(4.12) 



If q 

assertion of the proposition holds for cylinders). 



this follows from Lemma 4.16 (here we use the fact that we know that the 

> on D s . Note that for 



Now assume that q > 0. Then for all s < q we have 
all s < q we have D q C D s . Let R be the rectangle defined in Lemma 4.14 By Lemma 
4.15 (for s < q) we have 



sup 



p&D s C\D q C\R 



d_ 
dy 



sup 



p&D q C\R 



d_ 
dy 



-+0, 



as s — > q . This with the fact that for all s < q we have — 



> on D s implies that 



^ > on D q n R. But B 2 (D q ) cD g r\Rso Lemma 



dy 



We have shown (4.12). Note that by Lemma 



4.10 



4.12 



dy 

implies that 

dip, 



^ > on D„. 



dy 



we get 2p > o on W(D q 



It 



follows that > on D q . Now we will show that there exist a point po E B 2 (D q 



such that 



dy [Pi 



0) 



0. 



By (4. 12 ) we get that q < 1. By the definition of q we get that there exists a decreasing 
sequence {s n }^ 1 C [0, 1], s n — > q, such that for all n E IN we have (dip Sn /dy)(p n ) < 
for some p n E D Sn . 
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By Lemma 4.10 we may assume that p n G B 2 (D Sn ). Note that D Sn C D q . We also 



may assume (after taking a subsequence if necessary) that p n — > po G B 2 (D q ) 
also that p n <E R for any n = 0, 1, 2, 



■(Pn) - -^-(Po) 



<9y 



<9y 



< 



<9y 



< 



sup 

peD Sn nR 



dy 



. We have 








>-£<*> 


+ 




- dy M 




+ 


dy [Pn ' 


dy {Po) 



Note 



(4.13) 



By Lemma 4.15 the first expression in (4.13) tends to as n — > oo. The second 
expression in ( 4.13[ ) tends to as n — > oo because p n — ? ► pp, p n G D q D R (n = 0, 1, 2, . . .) 



and ^ is a Lipschitz function on D q n i? by Lemma 



Since ^ 



on D„ so ^ 



(Pn) 

(Po 



< it follows that ^(p ) 



4.14 



< 0. We know that 



dtl>q 

dy 



> and 



9lpq 

dr 



= 0. But po G B2(D q ), a contradiction with Lemma 

5. MONOTONICITY OF THE ODD EIGENFUNCTION 



4.8 



> 
□ 



In the previous section we proved monotonicity properties of for a special class 
of piecewise smooth domains Wi. In this section we pass to the limit to obtain the 



same result for class W. In the whole section we use notation (3.6). 



proof of Theorem 3.3 Note that by scaling it is sufficient to show the assertion of the 
theorem for D G T) such that r = 1. 
Let us consider the following inequality: 

V^Jfe) > V>( r i,yi) for an y 0i,2/i), (r 2 ,y 2 ) G D, r 2 > r u y 2 > yi- (5.1) 

Our first aim will be to justify (5.1) for all D G D such that r = 1. To show this 
inequality we will use the following scheme. If we have a family of sets {D s : s G [0, s }} 
such that all these sets belong to D(e, M, H, 1), W s = W(D S ), d W ( E ,M,H ,i)(W s , W ) -»■ 



as s — > and (5.1) holds for all D s where s G (0, So] then Lemma 3.16 guarantees that 
(5.1) holds also for Dq. We will show that (5.1) holds for all D G D such that ro — 1 



in a number of steps. 

Step 1. Note that Proposition 4.17 gives (5.1) for D G T> 1 . 

Step 2. Let T> 2 be the class of domains satisfying all conditions for the class T>i 
except that the condition (e) is replaced by 

(e') t —> r(t) is a continuous function on [Ti,T]. 



We will show that (|5.1| holds for D G £> 2 - Fix D G D 2 and let y , r x = T u T 2 , T 3 , 

min(ri,T — T 3 ,T 2 — Ti), = — yo + 1, so that 



T such as in Definition 



4.1 



Let £ 



-D G 1, .ff, 1), and let r(t), y(t), t G [0,T], be the parametrization of B(D). 

For s > let ft, 8 G C°°(R) be such that h s > 0, supp(/i s ) C (s,s), J R h s = 1. 
Put s = e/2. For s G (0,s ] let y s (t) = y(t) for t G [0,T], r s (t) = r(t) for t G 
[0,Ti+e/2] U [T -e/2,T] and r s (t) = f^ s r(t - u)h s (u) du for t G [Ti + e/2, T - e/2]. 
Let D = D and for s G (0,so], let D s be a domain defined like D but with (r(t),y(t)) 
replaced by (r s (t),y s (t)). Note that D s G D x for s G (0,s ] and D s G D(e/2,l,g,l) 
for s G [0, s ]. We have ||?/ s - y ||oo = and \\r s - r \ 



— )■ as s — )■ so Lemma 3.10 



ON HIGH SPOTS OF THE FUNDAMENTAL SLOSHING EIGENFUNCTIONS 



33 



gives that dyf( e /2,i,H,i)(W 8 ,Wo) — > as s — > 0, where W s = W(D S ). Hence Step 1 and 



Lemma 3716] give that (|5TT|) holds for D 
Step 3. 



D e D 2 . 



Let D 3 be the class of domains satisfying all conditions for the class D\ 
except that condition (e) is deleted and condition (c) is replaced by 

(c') t — > y(t), t y r(t) are strictly increasing continuous functions on [T 2 ,T 3 ]. 

It is clear that D 2 C D 3 , and since we can replace (r(t),y(t)) by {r (y -1 (y + t — Tx)) , y + 
t — Ti) for t G [T 2 ,T 3 ] (this is just a reparametrization), we see that in fact D 3 = D 2 . 

Step 4. Let D 4 be the class of domains satisfying all conditions for the class T>i 
except that condition (e) is deleted and condition (c) is replaced by 

(c") t — > y(t), t — > r(t) are nondecreasing continuous functions on [T 2 ,T 3 ]. 

Fix D G 2) 4 and let yo, r% — T\, T 2 , T 3 , T be such as in Definition Let e = 
min(rT,T - T 3 ), H = -y + 1, so that D G T>(s,l,H,l), and let r(t), y(t), t G [0,T], 
be the parametrization of B(D). 

For s G [0, 1) define r s (t) = r(t) and y s (t) = y(t) for t G [0, T 2 ] U [T 3 , T], and let 



r s (t) = (l-s)r(t) + sln + 2 . , /; 

2/.(t) = (l-sMt) + s(y +t-2i) 

for t G [T 2 ,T 3 ]. Let D s be the corresponding domain. Note that D s G T> 3 for s G (0, So] 
and D s G 2)(e, 1, if, 1) for s G [0, sq]. We have ||?/ s — yo||oo — > and ||r s — r ||oo — >■ as 
s — > so Lemma 3.10 gives that d wej ff,i)(W„ Wo) — )• as s — )• 0, where W s = W(Z) S ). 
Hence Step 2, Step 3 and Lemma 3.16 give that (5.1) holds for D = D G 2)4. 

Step 5. Let D5 be the class of domains satisfying all conditions for the class T>i 
except that conditions (b) and (e) are deleted, T 2 = 7\ and condition (c) is replaced by 

(c'") t — > y(t), t — > r(t) are nondecreasing continuous functions on [Ti,T 3 ]. 

Fix D G D5 and let yo, ri = T\ = T 2 , T 3 , T be such as in Definition 4.1 Let 
£ = min(ri,T - T 3 ), H = -y + 1, so that £> G T>(e,l,H, 1), and let r(t), y(t), 
t G [0, T], be the parametrization of B(D). 

Let s G [0, 1/2] and T 2 (s) = (1 - s)Tx + sT 3 . We define r a (t) = r(t) and = y(t) 
for t G [0, Ti] U [T 3 , T], r s (t) = n and y fl (t) = y + 1 - T x for t G [T 1; T 2 (s)], and let 



r s {t) 



T 3 -t 
1 - s 



y s {t) = {l-s)y[T 3 - 



n-t 



sy(T 3 



for t G [T 2 (s),T 3 ]. The corresponding domains D s are in D 4 for s G (0, 1/2], and D s G 
D(s, 1, if, 1) for s G [0, 1/2]. Since r(t) and are continuous, we have \\y a — j/o||oo — > 
and ||r s — ro||oo — > as s — > 0. Again, Lemma 3.10 and Lemma 3.16 together with 
Step 4 give that glj holds for D = D e T> 5 . 

Step 6. Let ©6 be the class of all domains from the class D such that tq = 1. Fix 
.D G D 6 and let e, M, if, T, y , r(t), be such as in Definition 3.1 Let s G [0,e/2]. 
We define r(t) = tandy(t) = ?/ forf6 [0,s], r(i) = 1 and y{t) = t-Tioxt G [T-s,T], 
and 



r s (i) = s + (1 — s)r 



t 



T -2s 



T 



y s (t) = -sH y 



T — 2s 



T 
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for t G [s,T — s]. Let D s correspond to r s (t) and y s (t) for s G [0, e/2\. Observe 
that D s G V(e,2M,H,l) for s G [0,e/2] and L> s G D 5 f or g G (0,e/2L_ Furthermore, 
ho -> and ||r s -r | 
D G D fi . 



Ws-yoi 



give (5.1 ) for .D 



— > as s — > 0, so Lemma 3.10, Lemma 3.16 and Step 5 



Step 7. Step 6 shows (5.1 ) for all D G 2) such that ro = 1. This implies the assertion 



of Theorem |3.3| but with weak inequalities instead of strict ones. By Lemma |4.6| this 
gives 1^ > on W and ^ > on W+. However 1^, are harmonic functions in W 

° ox — oy — T ox ' ay 

which implies §f > on W and §2 > on W+. Of course, this gives strict inequalities 

□ 



in Theorem 13.31 



proof of Theorem \l . 1\ (ii), (Hi). Elementary geometric considerations (we omit the de- 
tails here) give that any axisymmetric, convex, bounded liquid domain satisfying the 
John condition belongs to the class of domains D. Hence Theorem (ii) follows from 



Theorem 3.3 The inequalities in Theorem 1.1 (iii) follows from Theorem |1.1| (ii) by 



Lemma 4.6 The last sentence in Theorem 1.1 (iii) is an easy corollary of Theorem |1.1| 
(i). □ 



proof of Proposition 1.2. Let F — {(x, y, z) : x + z < R , y — 0} for some R > 0. Let 



9 G (7r/2, 7r) denote the angle between the free surface F(D) and the rigid wall B(D). 
It was proved in (formula (13.3)) that (see also |21j ) 

^(r, 0) = ij)(R, 0) (1 - u{R - r) cot 6 ) + o{R - r) as r ->■ R~ . 

Let us recall that we may assume that tp > on D. It is a standard result that ip is 
continuous on D and ^(0, 0) = 0. Of course tp attains maximum on F so ip must attain 
maximum on F(D). 

liip(R, 0) = then clearly ip(r, 0) does not attain its maximum at r = R. If ip(R, 0) > 
then ^(R, 0) = v cot 6oip(R, 0) < and therefore ip(R, 0) cannot be the maximum of 
</>(r,0). " 

It follows that ip{r, 0) attains its maximum inside the interval (0, R) and |^(r, 0) 
changes the sign on (0, R). Hence ip(r,y) has its maximum in the interior of F(D) and 
^:{r,y) changes the sign in D. □ 
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